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Abstract 

We analyse the periodicity of persistent currents in quantum spin Hall loops, partly covered with an s-wave 
superconductor, in the presence of a flux tube. Much like in normal (non-helical) metals, the periodicity of 
the single-particle spectrum goes from $o = hje to as the length of the superconductor is increased 

past the coherence length of the superconductor. We further analyze the periodicity of the persistent current, 
which is a many-body effect. Interestingly, time reversal symmetry and parity conservation can significantly 
change the period. We find a 2<i>o-periodic persistent current in two distinct regimes, where one corresponds 
to a Josephson junction and the other one to an Aharonov-Bohm setup. 
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1. Introduction 

In a seminal work from 1986, Biittiker and Klap- 
wijk discussed the flux sensitivity of a piecewise nor¬ 
mal and superconducting metal loop [T], as showed 
in Fig.[5a). The model they considered in order to 
describe such a system - a single electronic chan¬ 
nel with a linearized spectrum - is as simple as it 
gets, yet captures the characteristic features asso¬ 
ciated with persistent currents in mesoscopic loops. 
Indeed, in the Andreev approximation, and in the 
low-energy regime, microscopic details of the model 
hardly matter. As long as the length of the nor¬ 
mal region is much larger than the coherence length 
of the superconductor, the persistent current will 
have the familiar saw tooth shape, both in the nor¬ 
mal and superconducting regime. What changes 
though between the two regimes is the periodic¬ 
ity of the superconducting current with the applied 
flux. Following a simple calculation of the excita¬ 
tion spectrum, Biittiker and Klapwijk were able to 
show how as the length of the superconducting re¬ 
gion is progressively increased, the periodicity of 
the persistent current is halved, going from <I>o to 
$o/2, with $0 = hje the quantum of flux. The 
almost thirty years since the 1986 paper have seen 
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many exciting discoveries in the field of mesoscopic 
physics. One of them is the advent of topological in¬ 
sulators mill [5]. Of particular interest to us here 
is the case of quantum spin Hall insulators and their 
one-dimensional helical edge states, for several rea¬ 
sons. They first offer a new realization of ID Dirac 
physics, that goes beyond linearization of quadratic 
spectra at the Fermi points. Second, helicity, that 
is the locking of direction of spin with direction of 
motion, protects transport against time-reversal in¬ 
variant impurities. More precisely, single-particle 
elastic backscattering is forbidden by time-reversal 
symmetry, leading the community in a vast effort to 
better understand the effect of inelastic scattering 
in these systems IS 13 0 El m m na HI m HI HU . 
Third, the interplay of helicity and superconductiv¬ 
ity imposes a constraint on the parity of the num¬ 
ber of quasi-particles, or fermion parity (FP), in 
SNS junctions based on helical liquids, as the one 
depicted in Fig. nib)- This results in a so-called 
fractional Josephson effect [T3[ia[l9]i a ^q- 
periodic supercurrent. A few works have already 
discussed the physics of persistent currents in heli¬ 
cal rings, highlighting the effects of magnetic and 
non-magnetic impurities |20j . or the hybridization 
between edge states in narrow quantum spin Hall 
rings [2Ij . In the present paper, we revisit the 
analysis of Biittiker and Klapwijk in the context 
of quantum spin Hall insulators. In particular we 
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analyse the crossover between the normal persis¬ 
tent current and the supercurrent, as the length of 
the superconducting region is increased, and dis¬ 
cuss their periodicity. We argue that a constraint 
from time reversal symmetry doubles the period in 
the normal case, as compared to Ref. [T]. A similar 
effect is already known to occur in the supercon¬ 
ducting case, and was put forward by Zhang and 
Kane in Ref. [55]. 

The outline of the paper is as follows. In Sec. 
we start by discussing the model and the vari¬ 
ous symmetries that constrain the periodicity at a 
given edge, both in the normal and superconduct¬ 
ing regimes, and contrast the results with the case 
of a non-helical metal. Then, in Sec.[^ we comment 
on the total persistent current, when both edges are 
taken into account. Finally, in Sec.j^ we give some 
conclusions. 


2. Model and persistent cnrrent at a given 
edge 

2.1. Normal helical ring 

2.1.1. Single-particle spectrum 

We are interested in modeling the setup of 
Fig. 0b). To that end, we first restrict our analysis 
to a single edge, say the outer one. Before following 
in the footsteps of Ref. [T| and computing the exci¬ 
tation spectrum in the superconducting case, let us 
consider first the normal case. We model the outer 
edge by a ID segment of size L along the x direc¬ 
tion and impose periodic boundary conditions. The 
flux is included via the minimal substitution of the 
momentum operator p^ = —ihdx, 

Px -)-Px- qA (1) 

with q the charge of the particles and A = $/L a 
vector potential. In the following we take q = —e, 
for electrons. At energies much smaller than the 
bulk band gap, the helical states are well described 
by Dirac fermions, with the following single particle 
Hamiltonian 

n = vf{Px + eA)a3 , ( 2 ) 

where vp is the Fermi velocity and as is the usual 
Pauli matrix, acting on spin space. Right moving 
electrons have therefore spin up while left moving 
electrons have spin down (the situation is reversed 


at the inner edge). The single particle spectrum 
consists of two branches 

2tt / $ \ 

£±.n($) = ±hvFY V ^ j ’ 

with <I>o = hje. Corresponding eigenstates are of 
the form 

, (4) 

with x+ = (1,0)"^, x_ = (0,1)"^ and the mo¬ 
mentum kn = 27rn/L is quantized due to periodic 
boundary conditions. 

2.1.2. Excitations 

In the absence of the flux, $ = 0, we take the 
chemical potential to be at the Dirac point, that 
is, all states with negative energy are filled. Given 
the Hamiltonian of Eq. (0, nothing indicates that 
the spectrum is bounded from below, which would 
mean the ground state energy is infinite. Of course, 
in a real system, there is a natural cutoff scale, 
as the spectrum is bounded both from below and 
above by the bulk bands. However, at this stage 
one can let the cutoff go to infinity and renormal¬ 
ize the ground state energy without affecting the 
general physics which is given by the low energy 
excitations with respect to the Fermi sea. We then 
define the ground state |0)o such that 

c^|0)o = 0, n<0, (5) 

c±,n|0)o = 0, n>0. (6) 

and we impose that |0)o has zero energy. The index 
0 here serves as a reminder that the state is defined 
for $ = 0. There are two types of excitations on top 
of |0)o. One can create a particle in the conduction 
band, for instance „>o|0)oi or create a hole in the 
valence band, for instance c±_„<o|0)o. Importantly, 
particle and hole excitations are independent. We 
can also define a many-body Hamiltonian (still for 
$ = 0) as 

H = f dx : A{x)'H'i’{x) : , ( 7 ) 

Jo 

with 'I'(a;) = [^/i+(a:),'0-(a^)]^ a quantum field and 
where : ... : indicates normal-ordering, ensuring 
that indeed o(0|i7|0)o = 0. Among the many pos¬ 
sible excited states, some will play a particular role 
in the following. These are the states with a fi¬ 
nite number of particles, but no particle-hole exci¬ 
tations. They are obtained by filling positive en¬ 
ergy states or emptying negative energy states in 
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Figure 1: (a) Original setup as proposed in Ref. A piecewise normal and superconducting loop of size L comprising a single 
conducting channel is threaded by a magnetic flux Close to the Fermi points, there are four available modes: spin up or spin 
down right movers (blue solid lines) and spin up or spin down left movers (red solid lines). The normal region has length dj\j 
while the superconducting region has length dg- (b) A quantum spin Hall loop, partly covered by an s-wave superconductor 
(shaded region) and threaded by a flux $. Each edge hosts half the degrees of freedom available in setup (a), resulting in two 
helical channels physically separated by the insulating bulk. 


the following way. The state with Nj^± particles is 
defined by 


—1 

if Nj>0, (8) 

Nj-l 


|-^i)o = R Cj,-n\0)o ) 

n=0 

if Nj <0, (9) 

In the following we use the notation |Ab|_,iV_)o in 


order to refer to a state with right movers and 
N_ left movers. The energy of these many-body 
states is simply given by [53] 

£;(v+, = 0) = Y. +1) • 

( 10 ) 

This prompts us to introduce the so-called chiral 
current operators 

J±(a;) = :'i/'±(a:)V’±(a:) : , (11) 

as well as the particle current operator 

J{x) = vf ■ 'l'l(a;)(T 3 'l'(a:) : . (12) 

Next, we want to consider the effects of a finite flux 
$ 7 ^ 0 on the many-body states. There are several 
ways to look at the problem, some of which can 
lead to paradoxical claims, all of which are more or 
less related to the absence of a lower bound in the 
spectrum. The first observation is that the flux dis¬ 
turbs the vacuum. Indeed, at the level of the single¬ 
particle spectrum, all states move up and down 


when $ is varied. When one quantum of flux $o 
has been threaded, the single-particle spectrum has 
mapped onto itself, however the many-body state 
has changed, as an electron-hole pair has been cre¬ 
ated on top of the vacuum 


|0)c 


cYc.,o\0)o = \N+ = 1,N_ = -1)c 


(13) 

Using the definitions of Eqs. 0 and ( [T^ , one finds 
that this state carries a finite current 2vf/L. How¬ 
ever, we would like to compute the current and the 
energy for an arbitrary value of $ between 0 and $o- 
The problem is that these quantities are ill-defined, 
since our regularization scheme was introduced only 
for $ = 0. One way to circumvent this problem is 
to use perturbation theory. We assume that the 
field is turned on at some time to and compute the 
current at a later time t. We decompose H into 


H = Ho + e[t- to)Hi , (14) 


with 


Hq = vf [ dx , 

^0 


Hi — eA f dx J{x) 
Jo 


(15) 

(16) 


and treat Hi in linear response. From the Kubo 
formula we have 


{J±{x,t)) = {J±{x,t))o 
eA 

-I-— / dti dxi C^{x,t,xi,ti) , (17) 
« Jto Jo 
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with 

C^{x,t,xi,ti) = -ie{t-ti){[J±{x,t),J{xi,ti)])o . 

(18) 

Recognizing the famous anomalous commutators of 
bosonization 1^. 


{[Jaiq),Ja'i-q')])o = -Scr,a'5q,q' 

with Ja=±{q) = / 

Jo 


aqL 

27r 


(19) 


dx e ^’^^Ja{x), we finally ar¬ 


rive at 


{J±{x,t)) = - ± 


$ 

$n 


t > tn 


( 20 ) 


Eq. (201 shows clearly that the number of right 


(left) movers changes by -1-1 (—1) as one quantum 
of flux is threaded. This is a manifestation of the 
so-called chiral anomaly of (l-l-l)d quantum electro¬ 
dynamics [25] . In the context of condensed matter 
though, one is expected to provide a physical expla¬ 
nation to resolve the anomaly. One way is to think 
of an underlying ID lattice of which the Hamilto¬ 
nian of Eq. ([^ is only an approximation close to 
the Fermi points. The effect of the flux is then to 
push all states in one direction in the Brillouin zone, 
effectively converting one left mover into a right 
mover. However, for helical liquids, the more ap¬ 
propriate way of thinking is to involve the two edges 
of the setup, which exchange particles through the 
bulk as the flux is increased. A more detailed dis¬ 
cussion of the persistent current when both edges 
are involved can be found in Sec. [H We can also 
derive the correction to the energy coming from the 
anomalous current. One straightforward way is to 
realize that the current flows, in the presence of the 
induced electric field E{x,t) = —dtA{x,t). This 
results in an additional energy of 

Ew = [ dt' f dx {—e){J{x,t'))E{x,t') 

J tQ J 0 

^ 27r X ^ X- 27r / $ 

( 21 ) 

In turn, the ground state energy in the presence of 
the flux becomes 

Ptt 1 

E(iV+,iV_,$) = hvF-- Y, NiiN, + 1) 

27r X ^ X- 27r / $ 

-I- Hvf — {N+ - N_)— + hvF— — 

L $Q L y 


The energies of several states are represented in 
Fig. With Eq. ((22|), we can recover the usual 


expression for the persistent current at zero tem¬ 
perature 

m = If. (23) 

The approach we have used here is in our sense 
physically transparent and does not put too much 
emphasis on the unphysical infinite spectrum. In 
regular metals, with e.g a quadratic dispersion, the 
energy and the current can be computed directly 
by summing up over the finite number of occupied 
states below the Fermi energy. In the Dirac case, 
an analogous approach is available and consists in 
regularizing the infinite sum with a smooth cutoff, 
therefore mimicking the bandwidth. A detailed cal¬ 
culation, similar to the one leading to the Casimir 
effect, is nicely described in Ref. [20] . 

2.1.3. Persistent current 

Now, the question arises of what happens for 
large values of the flux. Increasing the flux requires 
more and more energy, therefore driving the sys¬ 
tem far away from the ground state. However, for 
some values of the flux (integer or half-integer mul¬ 
tiples of $o)) different many-body states become 
degenerate, as shown in Fig. and the system can 
in principle use these degeneracies to relax to the 
ground state. The consequence of this mechanism is 
a periodic persistent current. As was emphasized in 
Ref. [T] , the microscopic details of the model hardly 
play a role, and Eqs. (22) and (23) can, at the end 


( 22 ) 


of the day, in principle be used both for a normal 
and a helical channel. The crucial difference how¬ 
ever, will lie in determining which crossings are pro¬ 
tected by symmetry. Indeed, in the helical case, the 
degeneracies circled in Fig. [^ are protected by time 
reversal symmetry. To connect these states, one has 
to backscatter a right mover with spin up into a left 
mover with spin down, or vice versa. Barring the 
presence of magnetic impurities, matrix elements 
between these states are usually zero, and no gap 
opens. The other crossings concern states that can 
be connected by two-particle backscattering events 
which are not forbidden by time-reversal symmetry. 
As a consequence, the persistent current in the he¬ 
lical case has twice the period as in the non-helical 
case described in Ref. [T]. The system only returns 
to the ground state after threading $ = 2$o. In 
Fig-i we give several plots of the persistent cur¬ 
rent, both in the normal and helical case, for differ¬ 
ent starting points. The doubling of the periodicity 
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due to the topological protection against single par¬ 
ticle backscattering in helical liquids is the main de¬ 
parture from the results of Biittiker and Klapwijk. 
Our analysis does not however take into account 
effects of inelastic processes which can lead to tran¬ 
sitions between bands unnsziiisiin], and provide 
the system with a way to bypass the topological 
protection, relax to the ground state and restore 
the $0 periodicity. In the case where the length of 
the ring is small compared to the inelastic mean free 
path (a few /rm in typical quantum spin Hall insu¬ 
lators such as HgTe/CdTe [S] or InAs/GaSb [H]), 
the energy level broadening is small compared to 
the level spacing and the band structure given by 
Eq. is preserved |28) . 

If a measurement of the persistent current (in 
the normal state regime) were done on time scales 
larger than the corresponding inelastic scattering 
time (in the system under consideration) then the 
protection due to time reversal symmetry would be 
relaxed. The absence of this constraint would re¬ 
sult in a $0 periodicity of the persistent current. 
For InAs/GaSb quantum wells, with a Fermi veloc¬ 
ity of — 10^ m.s“^ [30], the inelastic scattering 
time is estimated to be of the order of 0.1ns. 



Figure 2: Plot of the energies of several many-body states 
£'(7V+, N—, $), with El = hvpJL. We have represented the 
lowest states with (a) W_|_ + N— = 0, (b) + N- = 1. 


2.2. Supercurrent in SNS junctions 

2.2.1. Excitation spectrum 

We now turn to the case where part of the topo¬ 
logical insulator is covered with an ordinary s-wave 




Figure 3: Plot of the persistent current for different starting 
points, in units of 1q = evpIL. (a) In the normal (non¬ 
helical) case there are two inequivalent starting points and 
the persistent current is $o periodic. The solid line corre¬ 
sponds to starting with N+ = N- = 0 at $ = 0 while the 
dashed line corresponds to starting with either Nj^ = 1 and 
N— = — 1 (positive current) or 7V+ = — 1 and N— = 1 (neg¬ 
ative current), (b) In the helical case, taking into account 
the constraint from time reversal symmetry, there are four 
inequivalent starting points and the persistent current is 2$o 
periodic. The blue solid line corresponds to = N— = 0 at 
$ = 0, the blue dashed line to either = 1 and N— = — 1 
or A^_|. = — 1 and 7V_ = 1, the red dashed-dotted line to 
7V+ = 1 and N- = 0 and finally the red dotted line to 
= 0 and N— = 1. 


superconductor which induces pairing by proxim¬ 
ity effect. We still focus for now on the outer ring 
of size L and assume that pairing is induced on a 
length ds- The remaining normal region has length 
dN- The crucial difference with the previous case 
is that, in a mean-field description of superconduc¬ 
tivity, electrons and holes are not independent ex¬ 
citations anymore. The system is well described by 
the following many-body Hamiltonian 

H=^j dx'ii'^{x)'HBdG{x)'i'{x) , (24) 


where now 4* = and 


T-t-Bdcix) 


fvpiPx + A{x)e 2 ieAa;/S, \ 

( A(x)e2“'^“'/^ -vf{Px - eA)a3) ’ 

(25) 
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is the Bogoliubov-de Gennes Hamiltonian, with 
A{x) = Ao9{x)9{ds — x), and 9{x) the Heaviside 
function. At the heart of Ref. [1], is the computa¬ 
tion of the spectrum for arbitrary values of c?Ar and 
ds- The spectrum has two branches given by the 
eigenvalues of Hbug in Eq. (25), and its periodic¬ 


ity with the flux interpolates between <i)o and $ 0 /2 
as the length of the superconducting region is in¬ 
creased. In the limit ds,diy ^ ^ with ^ = Hvf/Aq 
the coherence length of the superconductor, the 
single-particle spectrum reads 


e±,n(^) = 

dN 


1 , ^ ^ 
n -|— ± 2 — 
2 $0 


(26) 


Incidentally, such a spectrum is that of ID Dirac 
fermions on a ring of size 2dN, with twisted bound¬ 
ary conditions [3TJ[T9]. Building on Eq. (22), we 


could then write the following expression for the 
energy of a given many-body state 




i=± 


$ 

$0 


(27) 

Note that the expression is more compact than in 


Eq. (22) because in the absence of a flux the two 


consecutive Andreev reflections at the NS interfaces 
impose effective antiperiodic boundary conditions. 
However, due to the built-in particle-hole symmetry 
of the BdG Hamiltonian, ’’right” and ’’left” movers 
are not independent excitations anymore, and only 
one fermion number is needed to describe excita¬ 
tions in the junction ini [IS]. More precisely we 
have Nf = —N- = N. The final expression for the 
energy of a state with N quasi-particles is then 


E{N, $) = hvF 


2c? jv 


N + 2 


(28) 


We see that as the flux is increased by half a quan¬ 
tum of flux the number of quasi-particles in the 
junction is changed by one. This is another man¬ 
ifestation of the chiral anomaly, now for Majorana 
fermions, and is sometimes dubbed fermion parity 
anomaly in the context of topological superconduc¬ 
tors [TSj. The energies of several states are repre¬ 
sented in Eig.j^ A more precise, microscopic calcu¬ 
lation along the lines of Sec. 2.1.2 would lead to the 
same expression. Note that in the non-helical case, 
all many-body states are twice degenerate. Einally, 
the persistent current at zero temperature is still 
given by 

f) P 

im = ^. (29) 



Figure 4: Plot of the energies of several many-body states 
with Ed = hup/(2djv)- We have represented the 
lowest states with = 0, dil, ±2. 


2.2.2. Persistent current 
In the helical case, the crossings at <!> = ±$o/4 
between the states with N = Q and V = =f 1 are 
protected by fermion parity conservation dZlIISj. 
Crossing these values of the flux, the system can¬ 
not relax to the ground-state, and as a consequence, 
the persistent current is $0 periodic instead of $ 0/2 
periodic in the case of a non-helical junction. Eur- 
thermore, it has recently been put forward that 
time-reversal symmetry might double again the pe¬ 
riod [52]. Even though the analysis was originally 
done in the short junction limit, the simple argu¬ 
ments exposed in Sec. |2.1.3] together with the spec¬ 
trum of Eig. [^ can be used to understand the long 
junction regime. The idea is that the crossing at 
$ = 0 between the N = ±1 states is protected by 
time-reversal symmetry while the crossing between 
the N = ±2 states is not and could be opened by 
two-particle backscattering processes arising from 
Coulomb interactions. This would, once again, re¬ 
sult in a 2$o periodic persistent current. 


3. Total current including both edges 

We have seen that in quantum spin Hall sys¬ 
tems helicity and time-reversal symmetry impose 
new constraints on the many-body states, that can 
change the periodicity of the persistent current. In 
Table we summarize the situation. Of course, 
such an experiment with a flux threading would in¬ 
volve both edges of the quantum spin Hall insula¬ 
tor. The question then arises what is the behavior 
of the total persistent current given by the sum of 
the persistent currents at each edge. In the obvi- 
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System 

Spinful 

Helical 

Regime 

ds <i 

ds^i 


ds^ i 

Spectrum 

$0 

4>o/2 

$0 

$o/2 

Current 

$0 

$o/2 

a<I>o 



Table 1: Periodicity of the spectrum and of the persistent 
current in the spinful and helical case, in the normal regime 
(length of the superconductor dg much smaller than the co¬ 
herence length and in the superconducting regime (length 
of the superconductor dg much bigger than the coherence 
length ^). In the helical case, ck = 2 if the constraint of 
time reversal symmetry is taken into account, a = 1 if it 
is relaxed. Similarly, /3 = 2 if both the constraints from 
time reversal symmetry and fermion parity conservation are 
taken into account, = 1 if time reversal symmetry is re¬ 
laxed and finally /3 = 1/2 if fermion parity conservation is 
further relaxed. 

ous case where the width of the topological insu¬ 
lator becomes small, then time reversal symmetry 
would not protect the system against elastic single¬ 
particle backscattering (as an electron could change 
direction without flipping spin, simply by tunnel¬ 
ing from one edge to the other one) and the to¬ 
tal persistent current would recover the periodicity 
of its counterpart in the spinful metal. In reality, 
tunneling between edges can never be really sup¬ 
pressed and gaps will always open in the spectrum. 
What matters then is the time scale set by such 
gaps compared to the rate at which the flux is var¬ 
ied. Furthermore, even in the ideal limit of a de¬ 
coupling between edges, they are never truly inde¬ 
pendent. Indeed, they are both needed in order to 
explain physically the chiral anomaly, described in 
Sec. |2.1.2} In the normal regime, we have seen that 
as one increase the flux by one quantum of flux, a 
particle-hole excitation is created at the outer edge. 
Although no charge has been created, a left mover 
has been converted into a right mover, or, using our 
convention for helicity, a spin down has been con¬ 
verted into a spin up. This is only possible because, 
similarly to the Laughlin argument for the quantum 
Hall effect, the opposite process happens simultane¬ 
ously at the outer edge and the whole system acts 
as a spin pump [I71131]. 

4. Conclusions 

Following the ideas outlined by Biittiker and 
Klapwijk in Ref. [T] , we have revisited the question 
of persistent currents in quantum spin Hall rings, 
with or without proximity induced superconductiv¬ 
ity. We have found that contrary to a spinful metal. 


helicity together with time reversal symmetry puts 
additional constraints on the transitions between 
degenerate states and, as a result, can double the 
periodicity of the persistent current. 
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